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A phase damping reservoir composed by iV-bosons coupled to a system of interest through a cross- 
Kerr interaction is proposed and its effects on quantum superpositions are investigated. By means 
of analytical calculations we show that: i-) the reservoir induces a Gaussian decay of quantum coher- 
ences, and ii-) the inherent incommensurate character of the spectral distribution yields irreversibil- 
ity. A state-independent decoherence time and a master equation are both derived analytically. 
These results, which have been extended for the thermodynamic limit, show that nondissipative 
decoherence can be suitably contemplated within the EID approach. Finally, it is shown that the 
same mechanism yielding decoherence are also responsible for inducing dynamical disentanglement. 

PACS numbers: 03.65.Yz, 03.67.Mn, 03.65.Ta 



I. INTRODUCTION 



Environment-induced decoherence (EID) has been 
widely accepted as the mechanism responsible for the de- 
struction of quantum superpositions. A great amount of 
research in the last two decades has put in evidence - 
theoretically pj and experimentally Q - the potentiality 
of this effect in explaining challenging issues in quantum 
mechanics, as for instance the nonexistence of macro- 
scopic superpositions and the measurement process. The 
EID approach is based on the coupling of a system S 
with a many degrees of freedom environment, which is 
theoretically described in terms of microscopic quantum 
particles in equilibrium at temperature T. The coupling 
with the reservoir forces the system to evolve through a 
nonunitary, irreversible and, in general, dissipative dy- 
namics. 

One of the conceptual difficulties of the EID program, 
as it has been pointed out by some authors, concerns the 
irreversible energy transfer from the system to the reser- 
voir, which prohibits the application of this approach to 
conservative systems. Several alternative methods have 
been proposed to settle this problem. In general, these 
approaches defend the so termed self-induced decoherence 
(SID), whose cause is not attributed to external degrees 
of freedom (environment). Among them one may quote, 
for instance, proposals of a generalized Schrodinger equa- 
tion Q, assumptions about fluctuations of classical pa- 
rameters or internal degrees of freedom yielding decoher- 
ence 0, and nondissipative decoherence as a result of 
the quantum state averaging over a thermal distribution 
of velocities |{|. Also, references 0,0 contribute to this 
list. 

A recurrent result in the SID program is the reduction 
of the system dynamics to the following form of the well- 



known phase-destroying master equation 0, Q : 



dp 
Of 



(1) 



r being a specific parameter in each theory. By means of 
straightforward calculations in the eigenbasis {|-E n )} of 
H one may prove that the above master equation admits 
the following solution: 



p{t) = £ p nn ,e-^- E »'W h x 



cxp 
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where p nn > = (E n \p(0)\E' n ) and j = r. The purity, V = 
Tr[p 2 ], then reads 
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When r is constant, for example, one can observe an ex- 
ponential decay with time for the off-diagonal terms of p, 
implying in a decrease of the purity ('P(O) > 'P(oo)), even 
though the energy of the system remain unchanged, i.e., 
Tr[Hp(t)] = Tr[Hp(0)}. Clearly, the master equation 
JTJ suitably describes nondissipative phase decoherence. 
Recently, however, the relevance of SID as a physically 
meaningful process has been questioned [Tol ] . 

In this work we show that it is also possible to ob- 
tain phase decoherence within the EID scenario, with- 
out any additional assumption to the standard quantum 
formalism if we consider a diagonal interaction. Our ap- 
proach is based on a many-boson model in which the 
interaction between the components is given by the well- 
known cross-Kerr coupling a] ai Oj a,j . This assumption 
is grounded in the fact that this kind of interaction is 
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physically conceivable, i.e., two hybrid bosonic degrees 
of freedom such as ionic vibrational mode and the quan- 
tized field may interact effectively in this way as has been 
demonstrated recently Besides, several information 
transfer processes have been developed based on cross- 
Kerr interactions [l2^ in such hybrid systems, and sev- 
eral sources of decoherence without transition between 
the states of different quantum numbers - thus destroy- 
ing coherence without causing energy dissipation - have 
also been discussed in ^3|. These facts lead us to think 
that decoherence under this kind of phase channel will 
be a real problem to be considered. One can also do 
engineering of a bath of N identical ion traps each 
one coupled to the same field, such that they satisfy the 
following conditions of Ref. If: (i) the detuning of the 
atom with the field A is different from the multiples of the 
frequency v of the vibrational mode (A ^ kv, k integer) 
thus avoiding to excite the vibrational degree of freedom 
of each ion; (ii) in the large detuning limit (g <C A; g 
being the ion-field coupling) in order to achieve a dis- 
persive regime; and (iii) A <C v if the non-Lamb-Dickc 
regime is desirable. Although usually the energy dissipat- 
ing sources of decoherence may also be present (via typ- 
ical oscillator-oscillator position coupling), in the above 
engineered bath case, a nondissipative kind of decoher- 
ence is expected. By means of analytical calculations 
we show that the model indeed produces nondissipative 
decoherence, transforming quantum superpositions into 
mixed states through a Gaussian decay. 



In addition, we show that bipartite entanglement can 
be dynamically destroyed by such a reservoir. Dodd and 
Halliwell showed how open system dynamics destroys en- 
tanglement by means of the same mechanism that de- 
stroys interference 0|. Here, analytical formulas allow 
us to make the same identification rigorously. 



This paper is outlined as follows. In Sec. [H]the reser- 
voir is modeled system of N harmonic oscillators 
in a thermal equilibrium and the coupling with the sys- 
tem of interest is assumed to be given effectively by a 
cross-Kerr interaction. Exact results demonstrating the 
capability of the model in yielding decoherence for ar- 
bitrary quantum states are presented in Sec lIIII where: 
(A) analytical results demonstrating the Gaussian decay 
of quantum coherences and an estimate for the decoher- 
ence time are obtained for short times; (B) the equilib- 
rium regime is analyzed; (C) the thermodynamic limit is 
implemented; (D) two examples of decohering dynamics 
for some special quantum states, such as the cat state 
and a superposition of Fock states, are given, and (E) a 
phase-destroying master equation is derived in the regime 
of weak coupling. In Sec lIVI we demonstrate the dynam- 
ical disentanglement induced by the same mechanisms 
responsible for phase decoherence. Finally, in SecEl we 
relate our concluding remarks. 



II. THE MODEL 

In ref.jllj. Semiao and Barranco have derived an ef- 
fective Hamiltonian for a model of two-level trapped ion 
coupled to the quantized field inside a cavity. Under 
certain conditions, which allow for the rotating wave ap- 
proximation and assumptions on both a large detuning 
between the ion and the field (dispersive limit) and a 
small value for the Lamb-Dicke parameter, the authors 
find the following effective Hamiltonian 



H e f f = Q a\ ai + ui 7i\ 0,2 + A a\ a\ a3> 0,2 ■ 



(4) 



This Kerr-type coupling produces only phase interaction, 
not being able to introduce energy exchanges between the 
bosons. For this reason, this kind of phase coupling turns 
out to be a potential candidate for modeling nondissipa- 
tive reservoirs. 

Having been inspired by we consider the following 
many-bosons model 



H 



N 

E 

i=0 



huji rii 



N 

E 

i.j = 
i>3 



hgijUi rij 



(5) 



where denotes the number operator of the i-th bo- 
son. The coupling is clearly incapable to produce energy 
exchanges between any two degrees of freedom. 

Now we focus on the dynamics of a single boson, say 
the one with index "0" , in order to investigate the fea- 
tures of such an environment. For simplicity, we also 
assume that there is no interaction between the reservoir 
components. The Hamiltonian JSJ is then replaced by 



where 





= huj n Q , 




N 


Hr 


= ^ hu k n kl 




k=l 




N 


Hi 


= hg ok n n k 




k=l 



Initially, the joint quantum state is assumed to be 

P -I3H R 

p(0) = p s (0) ® - 

Tr e~P HR 



(0) 

(7a) 
(7b) 

(7c) 
(8) 



which indicates that all oscillators of the reservoir are 
in thermal equilibrium at a temperature (fc^/S) -1 . The 
system of interest is initially in an arbitrary pure state 
Ps(0). 

As it has been shown by Ford et al [lfj , the most gen- 
eral quantum Langevin equation, describing macroscopi- 
cally the dissipativc dynamics of a quantum system, can 
be realized with an independent oscillator (10) model. 
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Also, these authors have demonstrated that several tra- 
ditional models of reservoirs inducing dissipative decoher- 
ence can be mapped on an 10 model, which, accordingly, 
may be regarded as a fundamental model for the EID 
program. 

Our reservoir, however, may be regarded as being 
grounded on a type of generalized 10 model which im- 
plements nonlinear couplings. In fact, the original model 
for the ions, from which the effective cross-Kerr coupling 
given in Eq. Q has been derived in ref . 01 > contains the 
fundamental ingredients of an 10 model, including energy 
exchanges. 



III. PHASE DECOHERENCE 

The rather simple form of the coupling allows one to 
solve analytically the dynamics by means of straightfor- 
ward calculations, with no prior need for weak coupling 
or specific statistical assumptions. The exact result for 
the system density matrix in the Hq eigenbasis is written 
as 

oo 

Ps(t) =Y1 Pm, m 'C m , n At)e-^ m - m '^ ot \m)(m'\, (9) 

?n,m'— 

where p s = Tr R [p], p m>m ' = (m\p s (0)\m') and 



N 



Cm,m' {t) — J 



k=l 



\ — g-/3 htu k e -i (ra-m 1 ) g ok t 



(10) 



The purity of the system state, Vs = T^s[pS 2 }, is given 
by 



Vs(t) = ^2 \Pmm> 
m.m'—O 



|Cm,m' (£)| 



(ii) 



with 



/v 



\c m , m ,(t)\ 2 = n s i 



fc=i 



sin 2 


(m 


-m') g k t 




2 




sinh 2 


/3hui k 
2 





(12) 



Notice in Eq. I|ll|) that all the time and temperature de- 
pendences come from the function C TOjTO '. Also, in @ we 
see that C m m > is the only function depending on the cou- 
pling parameters between the system and the reservoir. 
Thus, it is natural to expect that this function carries all 
the dynamical informations about the correlations of the 
system with the environment and, therefore, about the 
decoherence of the quantum state. 

The analytical results lj§|)- l|12|) provide a rather inter- 
esting explanation for the well-known requirement of a 
spectral distribution for the coupling parameters gok- For 
resonant (gofc = g) or commensurate (<?ofc = h <?, h inte- 
ger) couplings, the time dependence in Ea. (|l(Jf) behaves 
with M integer. Then, for all instants sat- 



as e 



-iMg t 



for all m and m', indicating that the quantum state 
has just recovered its initial condition of purity, namely 
'Psitn) = Vs{0) = 1- This certainly is an unexpected re- 
currence for a supposedly irreversible reservoir. However, 
far apart from this very especial conditions, incommen- 
surate values of coupling inherent to any physically ac- 
ceptable spectral distribution will set irreversibility just 
by avoiding simultaneous r ccurrcnccs in Cm.m' • This is 
indeed ensured by the product in ljTT)|l . These statements 
arc demonstrated in Fig^ where a typical behavior for 
\Cm,m'\ 2 is numerically simulated. Notice that all ar- 
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FIG. 1: Log-linear plot for the function |C m , m ' | 2 as a function 
of (m — m')Qt for iV = 100 and f3hu} — 1. Two cases of 
couplings were simulated, namely: (a) a resonant spectrum, 
gou = n/V^V, and (b) a Gaussian spectral distribution, go k = 
Ui/Vn) exp [- - {k 2 ~ N \ o)2 ] • In (a) fl = 0.1 and in (b) = 

0.1125. In both cases, Q 2 = E fc 5ofe = 0.01 (see Eq.dJ) and 
ko — 50. Notice the long-time recurrences in (a). 



guments given in this section remain valid even for the 
simplest case of resonant harmonic frequencies (u>k = <jS). 
Since this assumption docs not change significantly the 
main features of the model, it will be considered hereafter 
for simplicity. 



A. The short-time regime 

In Fig.JQl one may clearly realize the presence of a 
short-time regime in which the collapse of the character- 
istic function occurs. Actually, in this section we shall 
show analytically that this is the very time scale of de- 
coherence. 

Let us consider the regime of short times, namely (m — 
fn') 9ok i<l. In this limit Ea. (|12f> can be replaced with 

( m _ TO ')2 t 2 /r 2 



isfying t n = 2i:n/g, with n € N, one has C„ 



1. 



\Cm,m' (t) I - 

where the decoherence time T£> is given by 
t d = 2Q- 1 sinh(/3ftw/2). 



(13) 



(14) 
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The frequency G, denned by Q 2 = J2k9ok' does n °t de- 
pend on the specific profile of the spectral distribution, 
since it contains the information only about the global 
sum, which is supposed to be normalized (the reason for 
that will be discussed latter). Indeed, in Eq. l|13|) (and in 
FigEJ one may see that the decay occurs even for a res- 
onant reservoir. Then, one may conclude that the only 
role played by the spectral distribution is to provide the 
character of irreversibility for the reservoir. 

Note, by Eo. (|13[) . that the decoherence time has been 
defined regardless of the difference m — m', even though 
the instant of decay is strongly dependent on this differ- 
ence. This choice attributes to td the role of an upper 
bound for surviving coherences, i.e., it corresponds to the 
time scale beyond which decoherence has indeed affected 
all elements of the matrix ps- 

Another important characteristic of our model is the 
Gaussian decay given by Eq. 1113 f) .which has already been 
found in other approaches [J, |5|- This behavior im- 
plies in a sudden decrease in the purity of the (arbi- 
trary) quantum state for times longer than td- In fact, 
well beyond this time scale the quantum state will have 
reached its equilibrium configuration. In this regime, 
since C m . m <{t) = \C m . m <{€)\ e**W, one may write 

C m , m '{t > t d ) = <5 m , m ', (15) 

indicating that the dynamics will transform 



As(0) = £ p m ,m'\m)(m'\ 



m,m'= 



with 



into 



with 



Ps(t >td) = J2 



V(t>T D )=J2p 



(16) 



(17) 



|m)(m| 



2 



(18) 



(19) 



Note that the structure of the coupling given by 
Eq. (|7c() automatically determines the pointer basis, i.e., 
the basis in which the density matrix becomes diagonal. 
This allows us to conjecture that this reservoir will fatally 
destroy the entanglement for a state like \n) \m) + \n') \m') 
(this will be demonstrated in Sec lIVjl . 

Another remarkable result derives from equations H17|) 
and ljT§|) : V(0) is always greater than V(t > td), indi- 
cating the occurrence of a base-independent decoherence. 
Furthermore, since no assumption was made about the 
pure initial state, the decohering process is also state- 
independent, td being the general upper bound for the 
time scale of such a process. 



This generality associated to td attests that this time 
scale is in fact an intrinsic feature of the reservoir model 
and not of the particular quantum state of the system. 
This is perhaps the most important difference regarding 
the usual time scales found in literature, which in general 
exhibit a dependence with d~ n , d being the distance in 
the configuration space separating the two components of 
the particular superposition state under investigation (for 
examples of such a state-dependent time scales see refer- 
ences f° r the f ree particle problem, and 0,0, for 
oscillators models) . In our model this dependence might 
be introduced by incorporating the "distance" m — m' in 
the definition given by Ea. (|14f> . However, we have opted 
to attribute a general character to tq. Concerning the 
dependence with other physical parameters, such as char- 
acteristic frequencies, temperature and h, our time scale 
agrees qualitatively with those of the references quoted 
above. For instance, all time scales predict a decrease 
in the decoherence time as the temperature increases. A 
precise quantitative agreement is not expected in princi- 
ple since we are working, in contrast to those approaches, 
with nonlinear interactions. 



B. The long-time regime 

Now we analyze the equilibrium situation (t — > oo). 
In this case, a direct inspection of Ea. i|12|) allows us to 
calculate a lower bound, which for any m and ml ', reads 



\a 



m,m' | lb 



sinh 2 (J3huj/2) 



1 + sinli (phu/2) 



N 



(20) 



It is easy to see, for N large but finite, that the lower 
bound given by Eq. I|20[l can range from arbitrarily small 
values, in the regime of high temperatures (flhuj <C 1), 
to values arbitrarily close to the unit, in the regime of low 
temperatures ((3Hlo 3> 1). The latter regime, for which 
\Cm.m'\LB corresponds to the one in which the dy- 
namics is practically unitary, without decoherence. On 
the other hand, the former, for which \C m<m i^ B — ► 0, 
corresponds to the regime in which the effects of deco- 
herence are accentuated. In FigEl these regimes can be 
identified clearly. 

The abrupt transition between these two limits (cen- 
tral region of Fig|2J) corresponds to a regime in which de- 
coherences still occurs, but only partially, since a finite 
value for \C m . m >\ LB indicates the existence of surviving 
off-diagonal terms in the quantum state. The critical 
temperature, i.e., the one beyond which the effects of 
decoherence become accentuated, may be estimated by 
imposing that \C mym '\ LB = 1/2 in Eq. (|2"U|) . The result is 
a function of N: 



~2 



arcsinh 



1 



V¥Jn - 1 



(21) 



For temperatures in the vicinity of {ksPcrit) 1 the quan- 
tum state must contain off-diagonal terms contributing 
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FIG. 2: Lower bound \C\ LB as a function of the dimensionless 
parameter fihw for N = 10 2 , 10 3 , and 10 4 (from left to right). 
The critical values /3 cr ithu), estimated from Eq. l21H for each 
N, are given by the vertical dashed lines. 



fco = 50, one may verify that Q — > 0.70617 as N — > oo. 
Also in this case, it may be verified that Q ^ 0.88912, 
this upper bound occurring for N = 2ko = 100. 

With these considerations we can reinterpret the re- 
sults given by Esa.l|13 |) -I|21 fl within the framework of a 
reasonable thermodynamic reservoir possessing the fol- 
lowing characteristics: 

i- ) Gaussian decay for the short-time regime; 

ii- ) Decoherence time given by td = 2 Q~ x sinh (/3 hui); 

oo 

iii- ) Characteristic frequency given by Q 2 = Y^, 5ofc> 

fe=i 

iv- ) Irreversibility induced by a quadratically normaliz- 

able spectral distribution (<?ofc); 

v- ) Lower bound |C| 2 ls = and T cr u = 0. 

Item v-), which can be directly obtained by taking N — > 
oo in Ea. (|20[) and Ea. (|21|) . tell us that decoherence will 
occur in the thermodynamic limit even for arbitrarily low 
temperatures. 



to an increase in the value of the purity in the regime of 
long times. 



D. Examples 



C. The thermodynamic limit 

In the previous sections we have discussed the gen- 
eral characteristic of our cross-Kerr reservoir for an ar- 
bitrary value of N. Now, we investigate the thermody- 
namic limit, which better describes the physical reality 
of a reservoir. 

Rigorously, the thermodynamic limit is implemented 
by taking N -> oo and V -> oo, with (N/V) -> 1. Be- 
fore considering this limit in our model, however, it would 
be necessary to replace the coupling <?ofc by its original 
physical form gofc / VV, which derives from the process of 
quantization of the electromagnetic field 0] . With this 
correction, the frequency Q defined in Ea. (|14f) would as- 
sume the form Q = yJ^Jj^jV. Then, in the resonance 

gofc = 12, for instance, we would have Q = Qy/N/V, 
which possesses a well defined thermodynamic limit 12. 

This simple example suggests that we can implement 
the thermodynamic limit in our formulas just by re- 
quiring that the coupling gofc be quadratically normal- 
ized, i.e., that Q 2 = X]fe5ofc ^ e finite. F° r instance, for 
the resonant case mentioned above, the assumption that 
gok = 12/a//V will automatically set Q = 12. In a less 
trivial case, in which the coupling is given by a Gaus- 
sian spectral distribution, the thermodynamic limit can 
be ensured just by assuming that 



n 



9Qk 



exp 



7r(fc-fc ) 2 
2N 2 



(22) 



In this case, one may verify numerically that the relation 
Q < 12 is always obeyed for any N. For instance, for 



Any initially pure state for the system of interest may 
be expanded as 



\i>o) = J2pH \n), 



(23) 



where |p(n)| 2 is the photon distribution defining the 
quantum state. The corresponding density operator is 
then written as ps(0) = |V'o)(V'o| an d its matrix elements 
in Fock basis reads p m ,m' = p(m) p* (m 1 ). Thus, Eq. Hll|) 
becomes 



oo oo 



V 8® = EE IpMI 2 W m ')l \C m . m ,(t)\ 2 . (24) 



m— ra'— 



Now we apply our model to investigate numerically the 
decoherence process for: (i) the cat state 



IV'o 



2 1 



= -2|«| 2 



(25) 



and (ii) the superposition of Fock states 

+ \m 2 ) 



IV'o) = 



V2 



(26) 



with mi 7^ fTT-2- For this Fock superposition, 
for which an analytical calculation yields Ps(t) = 
(l + \C mii7n2 (t)\ 2 ) /2, it is easy to verify the classical- 
ization just by applying formulas (|16|) and 1)18(1 . Indeed, 
the equilibrium value for the purity, Vs(oo) = 1/2, is 
directly associated to the statistical mixture ps (oo) = 
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(|mi)(mi| + |to2)(to2|) /2. For the cat state, on the other 
hand, at the equilibrium time we have a mixed state with 
< Vs(oo) < 1. In FigOHwe show numerical simulations 
for the purity of these states. The decoherence time is 
also plotted, attesting both the validity and universality 
of the estimate ()14|l . Notice that for SI 
t d ~ 1CT 8 s. 




FIG. 3: Log-linear plot for the purity of (i) a cat state (dashed 
curve), with a — 2.0, and (ii) the superposition (|1) + |2})/\/2 
(solid curve). For these simulations, N = 100 and j3Ku) = 
0.01. Here we have used the spectral distribution employed 
in Fig^(b). The vertical line stands for the decoherence time 
Q.T D = 1.13 x 10~ 2 , estimated from Eq.JTIJ. 



E. Master Equation 

The interaction Hamiltonian given by Ea. (|7c() can be 
written as 



Hi 



V s ® V R , 



(27) 



where we have defined: Vs — no and Vr = gok n%. By 
means of relatively simple calculations we then deduce 
a master equation for our model in the weak coupling 
approximation. The result, in the interaction picture, 
reads 



dt 



h 2 



%>Ps 



(28) 



where V s = (Vr) V s and £, = (Vj)/(U fl ) 2 . The reservoir 



expectations, (V^ 
by 

{Vr) 



= Tr 



Ar(0) Vr , are given exactly 



gok 



N 
k=l 



N 



(Vg> 



dot 



a /3huj 



fc=i 



(Vr) 2 . 



(29a) 
(29b) 



Using the spectral distribution given by Eq. (|22|l one can 
rewrite these formulas as 



(30a) 



where 



<) 2 



N 



e 0hu> + x 



A 1 (N) + (V R ) 2 , (30b) 



k=l 



VnN 2 



exp 



nN 2 



(31) 



One can verify numerically that A n (N) is a real con- 
stant obeying A n (N) ^ 1 for all N. Now we can imple- 
ment the thermodynamic limit just by taking N — ► oo, 
with some additional care. For the high-temperature 
limit (/3hoj — > 0), we see that these expectation val- 
ues explode. However, this regime is supposed not to 
be suitably described by the master equation, which pre- 
sumes weak coupling. For the low-temperature limit 
{fiYiLO — > oo) on the other hand, a solution for the for- 
mulas (|30l) is allowed to exist for N — > oo, in agreement 
with the analysis carried out in Sec llll (Jl Finally, we 
note that in this case we have (V^) — > (Vr) 2 , indicating 
the disappearance of some types of correlations. 

The structure of equation l|28|l may be promptly iden- 
tified to the phase-destroying master equation com- 
monly used in SID approaches. However, a rather impor- 
tant conceptual difference must be noted in the (nonuni- 
tary) diffusion term, which is in fact the one responsible 
for phase decoherence. In SID approach, H denotes the 
Hamiltonian of the system, whereas in EID approach the 
decoherence of the system is induced by its interaction 
with the reservoir. If on one hand these approaches share 
the same mathematical structure, on the other they can- 
not be conciliated physically. 

The dependence of Eq. (|2"8|l on t occurs due to the na- 
ture of the coupling employed to model the reservoir. 
In fact, this is expected for models of reservoir inducing 
Gaussian decay in the quantum coherences, as can be 
concluded from the formulas QJ-©- A natural extension 
of our model in order to contemplate the exponential de- 
cay commonly observed in several contexts may be tried 
by assuming an adequate temporal dependence for the 
coupling gofc- Research on this topic is in progress. 



IV. DISENTANGLEMENT INDUCED BY 
PHASE DECOHERENCE 

Consider now the case in which the system of interest 
is composed by two coupled oscillators. The Hamiltonian 
of the system is now considered to be 



H s = hu a n a + hu b n b + Hg ab n a n b . 



(32) 
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Adopting the same algebraic procedure of Sec lIIII we ob- 
tain the characteristic function 



JV 



Cm a ,m' a (t) - J] 1 _ e _ l3h „ k ^.[K-O^ + l-rK'^l' ' 



fc=l 



which is associated to the two-boson state 



(33) 



Ps(t) = £ p ma , m ' a C ma , m , a (t), 



-i{in a —m' a ) uj a t 



=0 i> b m b m b 



where 



x e-^-^^KmbXmXl, 



(TO a m 6 |p a6 (0)|TO^m' b ) 



(34) 



(35) 



Also, it is possible to show that 



\C„ 




(36) 



The analysis realized in the precedent sections concerning 
the short and long-time regimes can be also applied here 
with minimal changes. For the short-time limit we obtain 



\C ma , m , a W-e-t 

m b ,m' b 

where the decoherence time is now given by 



N r 



T ab — £ 



fe=l 



9ak (m Q -m' a ) + gbk (m b - m' b ) 



2 sinh(/3ftw/2) 



(37) 



(38) 



These expressions show how the decoherence process 
transforms 



PS 



(o) = E 



p ma , m ' a \m a m b ){rn' a m' b \ (39) 

= m b m b 
=0 



into 



Ps(t>T ab ) = y pm a ,™ a \m a m b ){m a m b \, (40) 

* * m h m h 



which has no off-diagonal terms in the Fock basis. Here 
again, additional care is required for the analysis of the 
critical temperature regime (for TV finite), in which |C| 2 
may assume finite values in the equilibrium, indicating 
the presence of surviving entanglement. Also, as it oc- 
curred in the precedent sections, for extremely low tem- 
peratures the degree of entanglement may remain prac- 
tically unchanged. 



Let us now consider the simplest entangled initial state 
for the system of two oscillators, namely 



(41) 



with n 7^ m. Since the interaction does not couple dif- 
ferent number states, the above vector plays the role of 
a two-qubit state coupled to the reservoir. In this case, 
we may follow the procedure adopted in ref . ppl ] , in which 
entanglement has been evaluated by means of two-qubits 
measures, as for instance, the negativity. This quantity 
has the following definition [2lL l22j : 



Af[p{t)] = 2max{0,-A„ es (i)}. 



(42) 



where \ n eg(t) is the negative eigenvalue of the partially 
transposed density matrix in the joint space £ ab - Follow- 
ing [20] we obtain analytically: 

W*) = -\^ c -z^ c z^)- ( 43 ) 

In the high-temperature and short-time regime we obtain 



^neg (^) 



* («-™) 2 t 2 /^ 



where 



T ab = 2Q- b 1 fanh(J3hw/2), 



(44) 



(45) 



with Q\ b = J2k(9ak + gbk) 2 /!- Equations (JHJ) and 102) 
show how entanglement disappears as the dynamics takes 
place, provided the temperature is higher than the criti- 
cal value given by Ea. (|21|l for the case of N finite. Also, 
Eg . (|44|) makes explicit that the disentanglement is in- 
duced by the same mechanism that produces decoher- 
ence, namely, the Gaussian decay in the characteristic 
function C{t). 



V. CONCLUDING REMARKS 

Inspired by a recent work |llj] we propose and theoret- 
ically investigate the properties of a reservoir composed 
by harmonic oscillators coupled to the system of interest 
by means of a cross-Kerr interaction. We have shown 
that this type of phase coupling is sufficient to guarantee 
decoherence without dissipation. 

Since the analysis have been carried out analytically in 
virtue of the simplicity of the model, the main features of 
our cross-Kerr reservoir could be clearly identified. For 
the regime of short times and finite temperatures, a state- 
independent decoherence time has been estimated and a 
Gaussian decay for quantum coherences has been demon- 
strated through the analysis of the characteristic function 

Cm,m' (0 ■ 

Concerning the coupling, we have concluded that the 
precise profile of the spectral distribution is relevant only 
in order to guarantee the irreversibility of the dynamics, 
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not being important to decide whether the decoherence 
may occur or not. In fact, the mechanism responsible for 
the decoherence has been shown to be the Gaussian de- 
cay (derived from the features of both the coupling and 
the thermal state), which is always present in the charac- 
teristic function for short times independently either on 
the regime of temperature or the number of oscillators in 
the reservoir. 

By analyzing the long-time dynamics for finite reser- 
voirs we have detected regimes of partial decoherence, 
in which the occurrence of surviving off-diagonal terms 
in the quantum state is expected. In this case, it was 
possible to identify a critical temperature beyond which 
decoherence is expected to be total. In the thermody- 
namic limit (N — > oo), however, the regime of partial 
decoherence no longer exists since T cr i t — ► 0. 

We have also derived a phase-destroying master equa- 
tion without statistical assumptions in the weak coupling 
approximation. The thermodynamic limit has been de- 
rived for the regime of low temperatures, indicating the 



occurrence of decohering dynamics under these condi- 
tions. These results suggest that the EID approach is 
able to produce most of the SID results with no need for 
additional postulates in quantum theory. 

Last but not least, we have verified a Gaussian decay 
of the entanglement associated to the a bipartite system 
coupled to our cross-Kcrr reservoir. The result, which 
has been obtained for arbitrary T and N , illustrates the 
role of disentanglement being induced by the same mech- 
anism responsible for the decoherence. 
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